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Abstract: The forests of northern Mexico and the southwestern United States have been subjected 
to warmer temperatures, persistent drought, and more intense and widespread wildfire. Tree- 
ring data from four conifer species native to these borderlands forests are compared with regional 
and large-scale precipitation and temperature data. These species include Abies durangensis, Pinus 
arizonica, Pinus cembroides, and Pseudotsuga menziesii. Twelve detrended and standardized ring-width 
chronologies are derived for these four species, all are cross-correlated during their common interval 
of 1903-2000 (r = 0.567 to 0.738, p < 0.01), and all load positively on the first principal component of 
radial growth, which alone represents 56% of the variance in the correlation matrix. Correlation with 
monthly precipitation and temperature data for the study area indicates that all four species respond 
primarily to precipitation during the cool season of autumn and winter, October—May (r = 0.71, 
p < 0.01, 1931-2000), and to temperature primarily during the late spring and early summer, January— 
July (r 0 —0.67, p < 0.01, 1931-2000), in spite of differences in phylogeny and microsite conditions. 
The instrumental climate data for the region indicate that warmer conditions during the January-July 
season most relevant to radial growth are beginning to exceed the warmest episode of the 20th century 
in both intensity and duration. The strong negative correlation between temperature and tree growth 
indicates that these four conifer species may be challenged by the warmer temperatures forecast in 
the coming decades for the borderlands region due to anthropogenic forcing. This information could 
constitute a baseline to analyze the impact of climate change in other regions of Mexico and the USA, 
where conifer species are of great ecological and socioeconomical importance. 


Keywords: climate warming; dendrochronology; tree-ring chronology network; climate response; 
correlation; principal component analysis 


1. Introduction 


Understanding climate variability and the impacts of climate change in natural pro- 
tected areas is important for the development of management strategies for conservation, 
restoration, and the sustainable use of natural resources. Natural areas protect biodiver- 
sity and may provide multiple environmental services of basic importance to society [1]. 
The sustainable management of natural areas is therefore potentially important to the 
communities living in or near by the protected preserves. 

A limitation to developing and implementing management strategies of natural pro- 
tected areas, in particular those located in the borderland area of the Mexican side, is 
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the lack of information about the dynamics of plant communities, annual radial growth 
and biomass production of dominant trees, scarcity of climatic information records, and 
the representativeness and quality of the available data. This information is needed to 
develop hydrological balances and to determine the water and other environmental ser- 
vices provided by these preserves and the impact that land-use changes may have on 
variables affecting the hydrological cycle and biodiversity [2]. Climatic scenarios for north- 
ern Mexico based on the Intergovernmental Panel on Climate Change (IPCC) models for 
the period 2000-2099 indicate an increase in aridity and more frequent droughts [3], with 
an estimated reduction in water availability of approximately 30% [4]. This reduction will 
affect those species with higher water demands, and it may increase tree mortality and 
disrupt ecosystem stability [5-7]. Warming temperatures have a direct effect on the rate 
of photosynthesis and respiration processes, increasing the rate of infestation, changing 
forest structure and biodiversity, and disturbing the carbon cycle and other ecological pro- 
cesses [8]. Understanding the effect of climate warming in forest species present in semiarid 
forest communities of the borderland area of Texas and northern Mexico is important to 
propose mitigation alternatives and to preserve the environmental services provided by 
those ecosystems, which may also contribute to ameliorate the effect of climate warming [9]. 

The Maderas del Carmen Natural Protected Area (MDCPA) located in northern 
Coahuila is adjacent to and east of Big Bend National Park (BBNP), Texas. The Cafion de 
Santa Elena (CSEPA) preserve is located immediately west of Big Bend NP along the Rio 
Bravo (Rio Grande) in Chihuahua, Mexico. These three protected areas share similar flora 
and fauna, and they constitute a bi-national corridor that protects regional biodiversity and 
provides habitat for endemic species native to the Chihuahuan desert [10]. 

Predicted changes in the amount and timing of precipitation along with further warm- 
ing may have a strong impact on plant species in the borderlands area of the Chihuahuan 
desert [11]. Lower precipitation and warmer temperatures have already caused bark-beetle 
infestations, wildfires, and tree mortality in the temperate forests of the southwestern 
United States during the last 50 years [12,13]. Understanding climate—vegetation relation- 
ships is important in the borderlands, where prolonged droughts may have a particularly 
detrimental impact on the native conifer species [7]. The objective of this study was to 
identify the precise seasonal precipitation and temperature signal present in tree-ring data 
of four conifer species native to the borderlands (i.e., Abies durangensis Martinez, Pinus 
arizonica Engelm., Pinus cembroides Zucc., and Pseudotsuga menziesii Mirb. Franco). We 
hypothesize that given the environmental and phylogenetic differences among the studied 
conifer species, they will have a different climatic response function and will be affected 
differently by climate warming. To what extent might projected seasonal climate change 
threaten these native forest species? 


2. Materials and Methods 
2.1. Location of the Study Area 


The natural protected areas BBNP, MDCPA, and CSEPA are located at the northern 
edge of the Sierra Madre Oriental (Eastern Sierra Madre) in western Texas, USA, and the 
state limits of Chihuahua and Coahuila, Mexico. The area is geologically diverse, and much 
of the main mountain mass is composed of Tertiary volcanic rocks covered by a layer of 
Cretaceous sedimentary rocks [10]. The topography of the area is complex dominated by 
rocky uplands separated by steep and deep canyons. This heterogeneous physiography 
creates varied environmental conditions that favor biodiversity [14]. 

The BBNP has an area of 324,000 ha in western Texas bordering Mexico along the Rio 
Grande/Rio Bravo; it has an elevation of 595 m along the river to 2575 m in the Chisos 
Mountain, and it is the southernmost range in the continental United States. The BBNP 
exhibits contrasted climatic conditions due to the range in elevation, which cause a wide 
variation in moisture availability and temperature. Annual precipitation ranges from 
150 mm in the arid alluvial plains to 430 mm at the Emory Peak at the Chisos Mountains. 
Summer temperatures reach near 40 °C, and subfreezing temperatures occasionally occur 
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in the wintertime. This variation in climate contributes to the extraordinary diversity in 
plant and wildlife habitats present in the park [15]. 

The MDCPA encompasses an area of 208,381 ha in the municipalities of Ocampo, 
Muzquiz, and Acufia, Coahuila, Mexico. Elevation varies from 500 to 3100 m and precipita- 
tion ranges from 100 to 400 mm at the MDCPA. Mean annual temperature varies from 18 
to 22 °C. These conditions favor the presence of different vegetation communities along 
the elevation gradient; thus, desert grasslands and shrublands dominate low elevations, 
and mixed-conifer forest dominate upper elevations. Dominant vegetation above 1600 m 
elevation is composed by a pifion-juniper community (1400 to 2700 m), pine—oak forests 
(1900 to 2600 m), and fir-cypress forest (1600 to 2700 m) with the presence of Psedutosuga 
menziesti Mirb., and Abies durangensis var. Coahuilenses [16]. 

The CSEPA has an area of 277,209 ha and is located at the international frontier at 
the northeastern margin of the Chihuahua state, bordering BBNP to the north. At an 
elevation range from 575 to 1840 m, the dominant vegetation communities are grasslands 
and shrublands at lower elevations and an oak—pinon community at higher elevations with 
the presence of Pinus cembroides, Quercus grisea, and Quercus arizonica [17]. 

Climate conditions in the CSEPA were obtained from weather stations located at the 
municipalities of Ojinaga and Manuel Benavides in the state of Chihuahua, Mexico, and 
Chisos in the BBNP, Texas, USA. Dominant climate varies according to the region, with 
a precipitation range from 287 to 571 mm at lower and upper elevations, respectively. 
Mean annual average temperature ranges from 19.5 to 26.0 °C, with an average minimum 
temperature of 10.0 °C in January and a mean maximum of 27.8 °C in June [15]. The 
hydrology of the region is very limited with watersheds yielding ephemeral streams 
draining toward the Conchos and Rio Bravo rivers. Dominant soils are shallow in depth 
and classified as lithosols, mollisols and entisols; they are composed of deep gravelly loam, 
which is well drained and non-calcareous with a low water-holding capacity [17]. 


2.2. Tree-Ring Chronology Development 


Twelve conifer tree-ring chronologies from the borderlands region were used in this 
analysis (Figure 1, Table 1). To develop the chronologies, the field sampling procedures in- 
volved obtaining two increment cores at breast height from selected trees with a Swedish in- 
crement borer and cross-sections from stumps or dead trees with a chainsaw in the MDCPA, 
CSEA, and neighboring watersheds. The tree species sampled were Abies duranguensis var. 
Coahuilenses, Pinus arizonica, Pinus cembroides, and. Pseudotsuga menziesii. The total ring- 
width data for BBNP, Texas, were downloaded from the ITRDB (https://www1.ncdc.noaa. 
gov/pub/data/paleo/treering / measurements /northamerica/usa/tx042.rwl (accessed on 
10 March 2022)). 
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Figure 1. The locations of tree-ring chronologies from the borderland are numbered and plotted 
(red triangles). The three large protected natural areas in the Big Bend region of the Rio Brave 
are also mapped along with the international frontier and the municipal divisions in Chihuahua 
and Coahuila. 


Table 1. Network of ring-width chronologies along the international border between Texas and the 
Mexican states of Chihuahua and Coahuila. 


Site Name Site Code Species Dated Series Dating 
1 Big Bend BIGTXPSME P. menziesii 101 1473-2014 
ABCTXABDU A. durangensis 17 1903-2004 
MDCMXPSME P. menziesii 30 1807-2004 
MDCMXPIAZ P. arizonica 24 1810-2001 
a MDCMXPICE P. cembroides 34 1827-2004 
MDCMXABDU A. durangensis 17 1797-2003 
Santa Elena ELEMXPICE P. cembroides 14 1725-2014 
Serrania del Burro SBUMXPIAZ P. arizonica 46 1798-2015 
Sierra Rica SRIMXPICE P. cembroides 53 1670-2015 
Pajaritos PAJMXPICE P. cembroides 13 1850-2014 
Santa Elena ELEMXPICE P. cembropides 14 1725-2014 
Sandillal de San Marcos SANMXPIAZ P. arizonica 18 1810-2014 
Namiquiapa NAMMXPICE P. cembroides 18 1627-2014 


1 https://www1.ncdc.noaa.gov/pub/data/paleo/treering / measurements /northamerica /usa/tx042.rwl (ac- 
cessed on 10 March 2022). 
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The polished cores and cross-sections were crossdated following standard dendrochrono- 
logical procedures [18]. Ring widths were measured with a Velmex measurement system 
to a precision of 0.001 mm, and the measured time series were screened for dating or mea- 
surement errors using correlation analyses in the computed program COFECHA [19,20]. 
The Signal Free method was used to compute the ring-width index chronologies [21,22]. 
However, the Signal Free approach requires large sample sizes for best results, and some of 
the borderlands datasets are based on relatively few trees. Consequently, we computed 
the chronologies used in this study as the robust mean value of the index series that were 
detrended and standardized based on an age-dependent spline function. The variance of 
the derived mean index chronologies was also then stabilized by the methods outlined by 
Meko et al. [23] and Stahle et al. [24]. A regional average ring-width index chronology was 
computed for each species and plotted for comparisons between the species. 


2.3. Correlation of the Tree-Ring Chronologies with Regional Precipitation and Temperature 


The 0.5° gridded monthly precipitation totals were obtained from the Global Precipita- 
tion Climatology Centre (GPCC) [25,26]. The 0.5° gridded monthly mean temperature data 
were obtained from the University of Delaware [27]. This database was chosen considering 
that the gridded climate data are derived from quality controlled station data, and that the 
full data reanalysis product (1901-2010) is well recommended for regional water balance 
studies, particularly considering that includes the period covered for most of the tree- 
ring series involved in this study. Other gridded precipitation datasets, e.g., ERA-Interim 
ECMWFE or CHELSEA (Climatologies at high resolution for the earth’s land surface area), 
are high-resolution databases, but they have extension limitations, providing information 
for the periods 1979-2019 and 1979-2013, respectively. On the other hand, Terraclimate has 
limited validation in data-sparse regions such as the Chihuahuan desert in northern Mexico. 

To compute the climate response of each tree-ring chronology, the gridded instru- 
mental precipitation data were first totaled, and the temperature data were averaged on 
a monthly basis for the study area (28.5-30.5° N, 102-107° W). The period of analysis 
extended from 1946 to 2000 because the instrumental precipitation and temperature data 
are limited before 1946, and some of the tree-ring chronologies end in the early 2000s 
when the core samples were collected. A significant correlation of at least (p < 0.05) for 
the common period of 1946-2000 was used to define the monthly climate signal in these 
tree-ring chronologies. 


2.4. Principal Component Analysis 


Principal component analysis (PCA) consists of a dimensionality-reduction method 
of a large data set used to minimize the number of variables, while preserving as much 
information as possible. To perform a PCA analysis, it is important to consider that the 
data set fulfills some statistical requirements; among them is that the correlation between 
variables be linear, the data set be free of outliers, and data be continuous. The tree-ring 
data set employed in this study fulfilled these requirements. Other dimensionally-reduction 
methods such as partial least square (PLS) or the principal component regression method 
(PCR) were not useful for this purpose, given that the PLS is particularly suited when 
the matrix of predictors has more variables than observations, whereas the PCR is used 
for estimating the unknown regression coefficients in a standard linear regression model, 
which was not the purpose of the analysis. 

In this study, PCA was used to investigate the common growth signal detrended and 
standardized chronologies of total ring width. The time series scores of the first principal 
component with the highest variance in the correlation matrix of ring-width chronologies 
was then correlated with gridded monthly precipitation and monthly mean temperature 
over North America to document the strength and spatial pattern of the monthly climate 
signal in conifer growth over the borderlands region. A time series perspective on the 
signal of climate variability and climate change was also derived by totaling the monthly 
precipitation for each grid point and then averaging it for the study area (28.5-30.5° N, 
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102-107° W) from 1901 to 2017. These regional precipitation data were then totaled for the 
October—May cool season and correlated with PC1 of the 12 conifer chronologies from the 
borderlands for the common period (1931-2000). The monthly temperature data were also 
averaged for the same area and for the seasonal period January-July and then correlated 
with the scores of PC1. 


3. Results 
3.1. The Borderlands Tree-Ring Chronologies 


The 12 total ring-width chronologies developed in the borderlands region are listed in 
Table 1. The shortest chronology is the A. durangensis series from Big Bend National Park, 
Texas, that extends from 1903 to 2004, and the longest is the P. menziesii chronology that 
extends from 1473 to 2014, which is also from Big Bend. Four chronologies were developed 
from the Maderas del Carmen, Coahuila, and the remaining chronologies were developed 
at various locations in Coahuila and Chihuahua (Table 1). All of these conifer chronologies 
were collected from arid forest sites that are highly sensitive to inter-annual and decadal 
climate variability [28,29]. 


3.2. Co-Variability of the Borderlands Tree-Ring Chronologies 


The results of the PCA based on the 12 borderland chronologies for the common 
period 1931-2000 are presented in Tables 2 and 3. These 12 chronologies are all well 
correlated, all load positively on PC1, and the PC1 scores are very strongly and significantly 
correlated with the soil moisture balance over northern Mexico and the southwestern US. 
PC1 explains 56% of the variance, while PC2 and PC3 explain 9.7 and 8.0% of the variance, 
respectively (Table 2). The loadings of each conifer chronology on PC1 indicate that the 
center of gravity for these 12 chronologies is the Maderas del Carmen Protected Area 
(Table 3), with the Douglas-fir from MDC having the highest loading. The loadings on PC2 
mostly discriminate the MDC chronologies from the others. 


Table 2. Principal components analysis of the 12 conifer chronologies from the Borderlands region of 
Texas and Mexico, 1931-2000. 


Comp Eigen Prop Var Cumvar 
PC1 6.717 0.56 0.56 
PC2 1.166 0.097 0.657 
PC3 0.963 0.08 0.737 
PC4 0.779 0.065 0.802 
PC5 0.598 0.05 0.852 
PC6 0.433 0.036 0.888 
PC7 0.391 0.033 0.921 
PC8 0.322 0.027 0.947 
PC9 0.193 0.016 0.963 
PC10 0.173 0.014 0.978 
PC11 0.158 0.013 0.991 


PC12 0.109 0.09 1.000 
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Table 3. Loadings of the 12 conifer chronologies on the first three principal components, 1931— 
2000. Chronology naming conventions: first three characters = site code, TX = Texas, MX = Mexico; 
ABDU = A. durangensis, PSME = P. menziesii, PICE = P. cembroides, PIAZ = P. arizonica. 


Chronology PC1 PC2 PC3 
ABCTXABDU —0.299 —0.212 0.242 
BIGTXPSME —0.271 0.372 0.171 
ELEMXPICE —0.277 0.232 —0.078 
MDCMXABDU —0.321 —0.198 0.357 
MDCMXPIAZ —0.284 0.349 0.324 
MDCMXPICE —0.321 —0.329 0.077 
MDCMXPSME —0.329 —0.296 0.164 
NAMMXPICE —0.284 0.349 —0.324 
PAJMXPICE —0.244 0.516 —0.181 
SANMXPIAZ —0.31 0.22 —0.575 
SBUMXPIAZ —0.228 —0.141 —0.575 
SRIMXPICE —0.28 —0.091 —0.502 


The correlations among the four separate species are also very strong, especially 
when based on the regional averages of two separate chronologies per species (Table 4). 
The regional ring-width chronologies for the four species are plotted in Figure 2. The 
inter-species correlations among these four species chronologies range from 0.56 to 0.73 
for the common period of 1903-2000 (Table 4). All four species recorded low growth 
during the severe drought of the 1950s, especially P. arizonica (Figure 2). Above average 
growth was recorded by all species during the 1980s (Figure 2). The strong loadings on 
PC1, the significant correlation among the species chronologies, and the synchronous 
sub-decadal tree growth anomalies in the 1950s and 1980s all suggest that the radial growth 
of the four conifer species in the borderlands region is driven largely by the same regional 
climate conditions. 


Table 4. The correlation matrix is presented for the regional average ring-width chronologies based 
on four conifer species for the common period of 1903-2000. 


1 Species ABDU PSME PIAZ PICE 
ABDU 1.0 
PSME 0.727 1.0 
PIAZ 0.607 0.738 1.0 
PICE 0.567 0.731 0.729 1.0 


1 All correlations are significant at p < 0.01. 
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Figure 2. Regional average ring-width index chronologies based on two collection sites per species are 
plotted for their common intervals. The Abies average is based on ABCTXABDU and MDCMXABDU 
(site codes defined in Table 1), which are correlated at r = 0.58 for 1905-2000 (a). The P. arizonica 
average is based on MDCMXPIAZ and SANMXPIAZ, correlated at r = 0.47 for 1860-2000 (b). The P. 
menziesii average is based on BIGTXPSME and MDCMXPSME, correlated at r = 0.42 from 1811-2000 
(c). The P. cembroides average is based on ELEMXPICE and NAMMXPICE, correlated at r = 0.62 for 
1725-2014 (d). The correlation matrix for these four species chronologies is presented in Table 4. 


3.3. Climatic Response of the Borderlands Tree-Ring Chronologies 

Each chronology was compared with monthly precipitation and temperature obser- 
vations averaged for the study area to carefully define the monthly correlations of each 
species with climate. The association between the chronologies and monthly precipitation 
and temperature from the prior August to the current August is plotted in Figure 3. Most 
of the chronologies responded positively and significantly to monthly precipitation of the 
previous September to the current February and from April to July. In terms of temperature, 
a negative response for January, March, April, June and July for most of the chronologies 
was found (Figure 3; Table 5). 
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Figure 3. The monthly climate response for the 12 conifer tree-ring chronologies from the borderlands 
region is plotted, based on the correlation between each chronology and the monthly precipitation 
totals (a) and mean temperature (b) for the study area from the prior August to the current August. 
The response of the A. durangensis chronologies is plotted in blue, P. menziesii in green, P. arizonica in 
black, and P. cembroides in red. The average monthly correlation for the 12 chronologies is also plotted 
(heavy black line) and with the significance thresholds (p < 0.05). The chronologies are listed in 
Table 1 and correlation coefficients are presented for all chronologies, climate variables, and months 
in Table 5. 


The precipitation response for most of the species is strongest in autumn-early winter 
and in spring. P. cembroides is most strongly correlated with precipitation in autumn and 
early winter preceding growth, and A. durangensis is most strongly correlated with spring 
precipitation (Figure 3, Table 5). None of the chronologies are significantly correlated 
with March precipitation totals for 1946-2000, and only four are significantly correlated 
with June precipitation. All 12 chronologies are negatively and significantly correlated 
with temperature in January and June, and the strongest seasonal response to temperature 
extends from May to July concurrent with radial growth (Figure 3, Table 5). Ironically, the 
strongest correlation between these 12 conifer chronologies and temperature was observed 
in June when the positive correlation with precipitation is weak and only marginally signif- 
icant for just four chronologies (Figure 3). The results suggest an interesting partitioning 
of the moisture and temperature response in these conifers of the borderland region of 
Texas and Mexico. The average moisture response for the 12 chronologies is strongest 
in autumn and early winter, but the temperature response is strongest in early summer 
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(June-July). With the possible exception of A. durangensis, these four species have a very 
similar climate response. 


Table 5. Correlation coefficients computed between each tree-ring chronology and monthly pre- 
cipitation and temperature data averaged for the study. All correlations were computed for the 
period 1946-2000 when the instrumental climate and the 12 tree-ring chronologies are most strongly 
replicated. Correlation values in bold are significant (p < 0.05). 


Precipitation 
Tree Ring Series 

pAug pSep pOct pNov pDic Jan Feb Mar Apr. May Jun Jul Aug 

ABCTXABDU 0.19 0.31 0.34 0.28 0.29 0.29 0.35 0.07 0.46 0.56 0.24 0.20 0.02 
BIGTXPSME 0.29 0.25 0.20 0.50 0.42 0.34 0.33 0.12 0.28 0.30 0.02 0.42 —0.09 

ELEMXPICE 0.31 0.18 0.36 0.46 0.22 0.25 —0.03 0.04 0.21 0.08 0.32 0.10 0.12 
MDCMXABDU 0.11 0.29 0.32 0.45 0.33 0.33 0.30 0.08 0.43 0.39 0.31 0.15 —0.02 

MDCMXPIAZ 0.42 0.46 0.34 0.26 0.13 0.29 0.32 0.14 0.26 0.28 0.03 0.24 0.25 

MDCMXPICE 0.26 0.35 0.41 0.33 0.24 0.32 0.30 0.14 0.32 0.39 0.17 0.28 0.09 
MDCMXPSME 0.22 0.32 0.31 0.44 0.30 0.25 0.34 0.05 0.31 0.39 0.29 0.17 —0.01 
NAMMXPICE 0.22 0.23 0.14 0.58 0.41 0.20 0.32 0.20 0.25 0.27 0.07 0.35  —0.07 

PAJMXPICE 0.26 0.21 0.32 0.51 0.25 0.01 0.17 —0.01 0.19 0.27 0.21 0.34 0.05 

SANMXPIAZ 0.32 0.35 0.34 0.52 0.43 0.22 0.21 —0.07 0.35 0.38 0.30 0.32 0.07 
SRUMXPIAZ 0.24 0.28 0.39 0.20 0.16 0.36 0.06 0.09 0.29 0.17 0.18 0.30 —0.01 

SRIMXPICE 0.50 0.49 0.45 0.47 0.17 0.18 0.31 0.02 0.27 0.16 0.19 0.44 0.10 

Temperature 

ABCTXABDU 0.14 0.28 0.28 —0.08 0.06 0.38 —0.06 0.36 0.31 0.40 0.50 —0.29 0.02 
BIGTXPSME 0.13 0.19 0.25 0.41 0.02 0.34 —0.22 0.40 0.33 0.44 —0.26 0.25 0.02 
ELEMXPICE 0.16 0.06 0.29 0.11 0.22 0.49 0.36 —0.10 0.03 0.07 0.47 —0.23 0.19 

MDCMXABDU 0.12 0.24 0.26 0.19 0.05 0.36 —0.05 0.37 0.30 0.36 0.56 0.28 0.00 
MDCMXPIAZ 0.29 0.29 0.53 —0.14 0.00 0.30 —0.12 0.27. —0.25 0.26 0.29 0.26 0.47 
MDCMXPICE 0.18 0.14 0.35 —0.07 0.00 0.39 —0.20 0.24 0.24 0.24 0.43 0.39 —0.25 
MDCMXPSME 0.16 0.24 0.34 —0.05 0.06 0.46 —0.14 0.36 0.31 0.28 0.52 0.33 —0.04 
NAMMXPICE 0.05 0.16 0.33 0.39 —0.10 0.33 —0.25 0.44 0.30 0.40 0.38 0.32 —0.10 
PAJMXPICE 0.14 0.12 0.14 0.01 0.24 0.39 —0.16 0.24 0.10 0.08 0.45 0.40 —0.06 
SANMXPIAZ 0.37 0.31 0.32 —0.10 0.21 0.41 —0.24 0.33 —0.18 0.37 0.62 0.39 —0.19 
SRUMXPIAZ 0.10 0.20 0.24 0.15 0.05  —037 0.03 —0.16 0.03 0.34 0.31 0.48 —0.21 
SRIMXPICE 0.47 0.43 0.39 —0.04 0.00 0.36 —0.33 0.13 0.22 0.33 0.57 0.55 0.36 


The time series scores for the first principal component of tree growth for the 
12 chronologies were correlated with gridded monthly precipitation and temperature 
over North America from the previous August to the current August to map the intensity 
and spatial pattern of the climate signal in the borderland tree-ring chronologies (Figures 4 
and 5). Note the strong and widespread correlation between PC1 of conifer tree growth 
from the borderlands region and gridded precipitation over Texas and northern Mexico 
beginning in October of the year prior to growth and extending most strongly through 
January (Figure 4). Note also the dipole pattern of positive and negative precipitation 
correlations most vividly evident in January (also seen in correlations between instrumental 
precipitation totals; Dettinger et al. [30]). The positive precipitation correlation over the 
borderlands extends into the summer concurrent with tree growth, but it is much weaker, 
especially in March (Figure 4). 
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Figure 4. PC1 of the Borderlands tree-ring chronologies correlated with monthly precipitation data, 
1946-2000. The time series scores for the first principal component of tree growth for the 12 detrended 
and standardized ring-width chronologies are correlated with gridded monthly precipitation totals 
over North America. The precipitation data were obtained from the GPCC with a grid spacing of 
0.5° [24], and the correlations were all based on the common period 1946-2000 when the tree ring 
and instrumental precipitation data are well replicated. 
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Figure 5. PC1 of the Borderlands tree-ring chronologies correlated with monthly temperature data, 
1946-2000. Same as Figure 3, except PC1 of tree growth is correlated with the 0.5 grid of monthly 
mean temperature from the University of Delaware [25]. 


The PC1 scores for the 12 chronologies are correlated with gridded monthly mean 
temperature in Figure 5. The tree growth correlation with monthly temperature is strong 
over northern Mexico and Texas in January, weakens substantially in February and April, 
and then is strongest in late spring and early summer (May-July). Note the very strong 
negative correlation between the PC1 scores and June temperature over the Rio Bravo region 
(Figure 5k). The correlation of the individual chronologies with June precipitation is quite 
weak for the study area centered over the Big Bend region (Figure 3, Table 5). However, the 
spatial correlations of the PC1 scores indicate that the region of highest correlations for June 
is shifted east of the study area to the lower Rio Bravo sector (Figure 4k), which is possibly 
related to the influx of moisture from the Gulf of Mexico early in the monsoon season. 

The dipole pattern of correlation with temperature from the borderlands into the 
Pacific Northwest is strongest from March through June. The correlations between the PC1 
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chronologies with precipitation and temperature of the prior April, May, June, and July are 
weak or even slightly negative for the study area (results not shown). This suggests that 
any carryover effects of precipitation or temperature anomalies from the monsoon season 
of the year prior to growth tend to be small or inconsistent. 

The instrumental precipitation and temperature data for the study area (28.5-30.5° N 
and 102-107° W) were totaled and averaged for the particular seasons that appear to be most 
important to conifer growth in the borderlands region (Figure 6) to gain some perspective 
on the potential for climate change to impact tree growth in these protected natural areas. 
The observed precipitation data are plotted for October-May and the temperature data are 
plotted for January-July for the full instrumental period of 1901-2017 (Figure 6). The PC1 
scores are plotted along with the seasonal precipitation and temperature data for 1931-2000, 
when the tree ring data are reasonably well replicated (Figure 6). Note that the instrumental 
precipitation and temperature data for the study are negatively correlated at r = —0.54 for 
these two partially overlapping seasons (Figure 6). 
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Figure 6. The instrumental precipitation data were totaled monthly for each grid point and then 
averaged for the study area (28.5-30.5° N and 102-107° W) for 1901-2017. These regional precipitation 
data were then summed for the October—May cool season and correlated with the PC1 scores for the 
12 conifer chronologies from the borderlands for the common period 1931-2000 ((a); PC1 in red). The 
correlation is r = 0.71 (p < 0.001). The monthly temperature data were also averaged for the study area 
and for the months of January-July (plotted for 1901-2017; (b)). The PC1 time series of tree growth 
was inverted and is plotted along with the regional temperature series. The correlation is r = —0.67 
(p < 0.001), prior to inverting the PC1 time series. Regional precipitation for October—-May has 
generally declined since the wet extremes in the 1980s and 1990s (a) while January-July temperatures 
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have increased over the region (b). Some of the tree-ring chronologies unfortunately end in 2000, 
but the PC1 time series does follow the low-frequency changes in both seasonal precipitation and 
temperature over the borderlands region. 


One implication of this strong association between the first principal component of 
tree growth and climate at the inter-annual and decadal scale is that conifer growth in 
the natural areas along the international frontier will suffer with the recurrence of severe 
drought episodes such as the 1930s, 1950s, 1980s, and 2010s (Figure 6a). However, the 
recent changes in January-July temperatures over the study area may be an even greater 
concern. In spite of large inter-annual variability, January to July mean temperatures have 
generally increased since the 1980s (Figure 6b) and will likely impact conifer growth if this 
seasonal warming continues. 


4. Discussion 


The development of tree-ring chronologies to analyze the climate variability of a basin 
or a larger region such as the borderland area of Texas, USA and Chihuahua and Coahuila, 
Mexico provides an excellent opportunity to understand the climatic response on conifer 
species growth and potential anthropogenic influence on climate variability. This advantage 
derives from the distribution of the chronologies in this large area and the involvement 
of conifer species with different ecological requirements that capture in greater detail the 
dominant climate variability and changes attributed to climate warming [31-33]. 

The natural protected areas BBNP, MDCPA, and CSEPA are part of the Chihuahuan 
desert and are dominated by particular environmental conditions that give place to special 
microclimates, hydrological behavior, and biodiversity [10,34]. However, overall, the 
climate on this extended region seems to be homogeneous due in part to the dominance 
of semiarid conditions and the common impact of ocean—atmospheric modes [9,24,35-38]. 
Climate warming, however, may put at risk the great biodiversity of this region, particularly 
when there is consensus about the potential impact of climate change on the stability of 
semiarid ecosystems present in the southwestern United States and the Chihuahuan desert 
due to the warm-season vapor pressure deficit largely controlled by temperature and 
cold-season precipitation [39,40]. 

The footprint of climate variability, however, is not similar in all semiarid or in more 
humid ecosystems around the world; recent research shows that temporal changes in 
drought conditions and physiological responses differ significantly among forest ecosys- 
tems, species composition (broadleaved, conifers), position along ecological gradients, and 
dominant climatic conditions [41-43]. It has also been confirmed that the annual growth of 
some conifer species is correlated with the meteorological variations, extent, and frequency 
of droughts around the world [44]; however, there is consensus that tree mortality is in- 
creasing due to drought and heat in forest ecosystems worldwide [5]. Previous studies 
have attributed an increase in fire frequency, reduction in forest productivity, increase in 
insect outbreaks, and change in species composition, among other negative impacts on the 
ecosystem, to the influence of climate change [12,39]. 

The four conifer species investigated represent a range of ecological conditions. Pinus 
cembroides thrives in semiarid climates (300 to 700 mm) with relatively high temperatures 
and evaporation rates, in soils nearly leveled or steep slopes with shallow soils on south 
and east-facing slopes [45]. On the other hand, Abies durangensis var. Coahuilenses and 
Pesudotusga menzieii are part of a mixed conifer forest in MDCPA found at higher elevations 
(over 2300 m) on north-facing slopes and humid canyons, whereas Pinus arizonica is part of 
a pine forest community growing at similar elevations but in drier sites [10]. Nevertheless, 
despite these ecological differences, they responded similarly to the dominant climatic 
conditions of the region, with significant positive response to the seasonal cool season of 
autumn and winter (October—May) precipitation and to temperature during the late spring 
and early summer seasons (January-July). In spite of these environmental gradients and 
phylogentic differences among Abies, Pinus, and Pseudotsuga, the climate response is quite 
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similar in this region. These results may not differ from previous findings obtained with 
particular species present in similar ecological environments, but they are relevant from 
the climatic response perspective considering that the all the species involved in this study 
showed a common climate signal, where 56% of the variance was explained by the first 
principal component. 

The impact of climate warming in semiarid ecosystems of the southwestern United 
States and northern Mexico may have important implications for future forest dynamics, 
plant cover, and environmental services they may provide [9,11,44,46,47]. 

The October-May accumulated seasonal precipitation to which the borderland 
chronologies responded has shown a significant decline since the wet extremes in the 
1980s and 1990s, while January-July temperatures have increased over the region (Figure 5). 
If this temperature-raising trend continue to escalate, the vapor pressure deficit largely 
controlled by winter-early summer temperature and cold season precipitation will exacer- 
bate forest drought stress, causing species mortality and triggering other negative impacts 
on semiarid ecosystems [7,13,39,48]. This negative impact has already been forecasted for 
the southwestern United forests where forest species dieback [5], insect outbreaks [49], 
increases in tree mortality [50], and fragmentation of ecosystems have been underway [51]. 

In the 20th century, the 1950s drought with a low October—May precipitation and 
high January-July temperature for eight consecutive years (1950-1957) affected extensive 
areas of the United Stated and most of the republic of Mexico [52]. This drought was 
one of the most severe in the last 500 years and could be considered a useful analog for 
the environmental consequences of higher temperatures and lower precipitation that are 
affecting forest ecosystems in semiarid regions such as the borderland area [39]. 

The 1950s drought in Mexico had serious socioeconomic impacts, i.e., failures in crop 
yields, lack of grains, famine, and livestock deaths, among others [53-55], but there is 
limited information on the impacts that this drought produced on the ecosystems. In the 
case of the United States, the 1950s drought in near three-fourths of the conterminous 
United States produced water shortages and streamflow reduction for different productive 
activities [56]. The years 1951, 1953, and 1956 had annual precipitation <225 mm, while the 
wet 1980s and 1990s had over 300 mm in the state of New Mexico with a similar effect for a 
great part of the semiarid regions [57]. 

Recent studies based on General Circulation Models (GCMs) project an increased risk 
for megadrought development over the southwest due to greenhouse gas forced reductions 
in precipitation and increased temperature and evaporative demand [58]. These predicted 
climate changes in the borderlands region are expected to aggravate harmful insect out- 
breaks [49], cause greater tree mortality [5], increase fire frequency and severity [47-49], 
change species composition [5], and lead to ecosystem fragmentation [44]. These projected 
changes threaten the stability of the conifer ecosystems in the protected natural areas of the 
borderlands region and the environmental services they provide [40]. 

As climate change continues to increase in the borderland area as projected by climate 
models, the mean forest drought stress by 2050 will exceed that of the most severe droughts 
on the last millennium [39]; thus, the species distribution, biodiversity, and dominant 
ecosystems present in this region will be altered. An understanding of the trends in these 
changes can allow protected area managers to implement more effective climate-adaptation 
strategies to conserve some of the most sensitive species to climate warming. 


5. Conclusions 


Four conifer species dominant in natural protected preserves of the borderland area 
with different habitat types responded similarly to the regional average total precipitation 
October-May and to the January-July average monthly temperature. Regional precipitation 
for October—May has generally declined since the wet extremes in the 1980s and 1990s, 
while January-July temperatures have increased over the region. The simple fact that 
these entire chronologies load positively at about the same magnitude on PC1 is telling 
us that they are very similar in their climatic response. This suggests that tree growth in 
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the region is mainly driven by moisture, but the increased trend in temperature observed 
at the regional level may exacerbate the 1950s drought and put at risk the ecosystems of 
these borderland-protected areas. These results do not contradict climate response results 
previously reported for some of these species from the greater southwestern sector of North 
America, but they do have important implications for future forest dynamics, plant cover, 
and the environmental services they may provide. 
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